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ABSTRACT 

We examine the faint-end slope of the rest-frame V-band luminosity function 
(LF), with respect to galaxy spectral type, of field galaxies with redshift z < 
0.5, using a sample of 80,820 galaxies with photometric redshifts in the Cosmic 
Evolution Survey (COSMOS) field. For all galaxy spectral types combined, the 
LF slope ranges from —1.24 to —1.12, from the lowest redshift bin to the highest. 
In the lowest redshift bin (0.02 < z < 0.1), where the magnitude limit is My < 
—13, the slope ranges from a ~ —1.1 for galaxies with early-type spectral energy 
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distributions (SEDs), to a ~ —1.9 for galaxies with low-extinction starburst 
SEDs. In each galaxy SED category (Ell, Sbc, Scd/Irr, and starburst), the faint- 
end slopes grow shallower with increasing redshift; in the highest redshift bin 
(0.4 < z < 0.5), a ~ —0.5 and —1.3 for early-types and starbursts respectively. 
The steepness of a at lower redshift could be qualitatively explained by large 
numbers of faint dwarf galaxies, perhaps of low surface brightness, which are not 
detected at higher redshifts. 

Subject headings: galaxies: luminosity function, mass function - - galaxies: 
evolution — galaxies: dwarf — galaxies: fundamental parameters — cosmology: 
observations — surveys 



Introduction 



The luminosity function (LF) of galaxies varies substantially with respect to many key 
physical parameters such as galaxy morphology, environment, color, star formation rate, 
surface brightness, and redshift. These many differences serve as powerful diagnostics of the 
broad tapestry of galaxy evolution. 

The most important requirement for the accurate derivation of galaxy LFs are large, 
complete samples of galaxies with reli able photometry and redshift determinations. The 
Cosmic Evolution Survey (COSMOS; IScoville et al.l l2007al ) contains not only the largest 
contiguous area of the sky yet observed with the Hubble S pace Telescope, bu t also deep 
multiwavelength imaging across the entire 2 □''COSMOS field JCapak et aDl2007h . COSMOS 
thus affords deep, homogeneous photometric and photometric redshift coverage for a sample 
of some 10 6 galaxies - complementing well at higher redshifts the largest galaxy surveys of 
the relatively nearby universe with which galax y LFs have been derived, such as the Two- 
Degree Field Galaxy Redshift Survey f2dFGRS: IColless et alihoOll: ICroton et al.l l2005h and 



the Sloan Digital Sky Survey (SDSS; iBlanton et al 



2003 



Abazajian et al.ll2004l ). 



Comprehensive analyses of the LF characteristics of the entire COSMOS galaxy sam- 
ple will ulti mately be forthc oming;, with the completion of the spectroscopic portion of 
the survey (ILilly et al.l 120071 ) and the continued addition and refine ment of the photo 



metric and photomet ric redshift measurements in multiple bandpasses (jCapak et al. 



2007 



Mobasher et al.ll2007l ). Already, though, it is feasible to address important scient ific questions 



abou t galaxy LFs with the current optical and near-infrared mult iband data (jCapak et al. 



20071 : iTaniguchi et al.l 120071 : lAjiki et al.ll2007l : IScarlata et al. 



2007) in conjun ction with the 



Hubble Space Telescope 7 8 i4 broad-band ACS imaging (jScoville et ah 2007b ). 
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One particular component of the field galaxy LF, the faint end (Mab > —18) of the LF 
at low to intermediate (0 < z < 0.5) redshifts, is uniquely well-suited for analysis with the 
COSMOS data. This population of galaxies lies in an apparent magnitude range somewhat 
beyond the 2dFGRS and SDSS survey limits; they are too faint for spectroscopic observa- 
tions, with redshifts too low to be included in surveys optimized for high-redshift galaxies. 
Yet, exploring these galaxies' contribution to the overall LF is critical for understanding 
galaxy evolution during the latter half of cosmic history. 

Recent comprehensive studies of field galaxy LF evolution have focused on the luminosity 
evolution of the bright end (Mab ^ —20). There, some consen s us appears to be grad ually 
emerging about the ex tent of that evolution. iDahlen et al.l (120051 ) . IWillmer et al.l (120061 ). and 
Scarlata et al. I (120071 ) all find a brightening of M B * by ~ 1 magnitude between < z < 1. 
In shorter-wavelength bandpasses, the luminosity evolution is more pronounced, while at 
longer wavelengths it appears to be weaker and p ossibly even in the negative sense, dimming 
in the near-infrared J-band from z ~ 0.4 to 0.9 (IDahlen et al.ll2005l ). 



The evolution of the faint-end slope, however , remains highly uncertain. At low red - 
shifts, results from the Sloan Dig ital Sky Survey (IBlanton et al.l 120051 ; iBaldry et al.ll2005l ). 
the 2dF Galaxy Reds hift Survey (Croton et al. 20051 ) , and other data sets such as that from 
the Century Survey (IBrown et al.l l200ll ) and the GALEX mission (IBudavari et al.l 120051 ) 
roughly agree, for example, on a moderate slope of a ~ —1.1. Beyond redshifts of a few 
tenths, the faint-end slope becomes very difficult to address, mainly because the number of 
low-luminosity galaxies detected in galaxy surveys decreases dramatically with increasing 
redshift. Despite a number of efforts to measure evolution in the faint- end slope at redshifts 
less than z ~ 1 (jWolf et al.ll2003l ; lllbert et al.ll2005l ; IZucca et al.ll2006l ). very little is known 
for galaxies fainter than M ~ —18. This is in part because the relationship between the 
bright-end and faint-end characteristics of galaxy LFs is not straightforward, o ften resulting 



i n a t radeof f between the p r ecisio n LF evolution measurements at the two ends. IBaldry et al. 



( 120051 ) and IWillmer et al.l (120061 ). for example, each select fixed LF faint-end slopes based 
on galaxies brighter than M ~ —18, and use that constraint throughout their bright-end 
evolution measurements. 

Whether or not the faint-end slope evolves with redshift, however, it is clear that its 
steepness varies widely for galaxies of different morphological and spectral types, indicating 
substantial differences in the evolutionary histories of galaxies. Broadly speaking, irregular 
galaxies, "blue" galaxies, and strongly star-forming galaxies - three significantly overlapping 
galaxy sub-populations - evolve more strongly than galaxies of other types. Such galaxies 
are characterized by very steep faint-end LF slopes and subs tantial evolution in luminosity 



and/or number density at even moderate redshifts of z < 0.5 (IMarzke et al.lll994l ; iLilly et al. 
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199a Ellis et al.lll996l : iLiu et al.lll998l : iBromlev et al.lll998l : iLin et al.lll999l ). More recent 



Perez- Gonzalez et al. 


2005fc 


Dahlen et al. 


2005) 



that composite parameteriz ations may more accurate l y reflect the shape of the L F than the 



usual single-function ones (Ide Lapparent et al.l 120031 . 12004 ; iBlanton et al.l 120051 ); and that 
faintward of M ~ —18, the powe r- law slope of the LF may differ quantitatively from the 



slope brightward of that threshold (IMadgwick et al.ll2002l ; iNorberg et al.ll2002l ; IBlanton et al. 



2005h . This could arise, for example, due to large nu mbers of low-luminosity " blue" galaxies 
(IWolf et al.ll2003l ) or low surface brightness galaxies (jlmpey et al.lll996t IBlanton et al.ll2005l ) 
which may have previously evaded detection. 

In this paper, we present measurements of the faint-end slopes of the rest-frame V- 
band luminosity functions of galaxies in the COSMOS survey at < z < 0.5, focusing in 
particular on the change in that slope as a function of redshift and galaxy spectral type. The 
depth and breadth of the COSMOS multiband photometry allows for reliable identifications 
of galaxy redshift and spectral type, with robust redshift error estimates for each galaxy, 
to a limit of tuab ~ 25 in the optical passbands. Even so, the relatively large and varying 
redshift uncertainties o f photometric redshifts can present substantial quantitative challenges 



and s ystematic biases (ISubbaRao et al.lll996t ILiu et al.lll998t IChen et al.ll2003l ; iDahlen et al. 



20051 ). We use Monte Carlo simulations to characterize these biases, and to recover the faint- 
end slopes of galaxy-type specific LFs. 

This work represents an initial study of the general properties of these type-specific 
faint-end LF slopes, to provide quantitative comparisons with the results of other large field 
galaxy surveys. Studies of the bright-end evolution o f the galaxy LF from the COSMOS 
survey are given elsewhere (e.g., IScarlata et al. 1 120071 ); and a more detailed breakdown of 
the < z < 0.5 galaxy population by redshift, galaxy spectral type, and galaxy morphology 
will be presented in a future paper (Liu et al., in preparation). 



Throughout this paper, we adopt a flat cosmology with Q\ 
= 70 km s- 1 Mpc" 1 . 



0.7, n r 



0.3, and 



2. Galaxy Sample and Photometric Redshifts 



I n our analysis, we use a compilation of the COSMOS optical/near-infrared data (ICapak et al. 



20071 ; iMobasher et al.l 120071 ) . which includes observations with HST Advanced Camera for 
Surveys (hu), the Subaru Telescope (B, V, r', i', z', NB816), the Canada- France-Hawaii Tele- 
scope (u*, £*), and the 4-meter KPNO Mayall and CTIO Blanco Telescopes (K s ), as well as 
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supplementary data from the Sloan Digital Sky Survey. The data from the different tele- 
scopes were all match ed to a common pixel sc ale and smoothed to the same point source 
function. SExtractor (IBertin fc Arnouts Ill996l ) was then used in dual mode to generate a 
photometric catalog, selected using the Subaru %' and CFHT i* images. The limiting 3a 
AB magnitude is i' = 26.03. A det ailed descript i on of the imaging data, photometry, and 
photometric calibration are given in ICapak et al.l (120071 ) . 



Photo metric redshifts for in dividual galaxies were computed using the methods de- 
scribed in iMobasher et al.l (120071 ). Six basic gal axy spectral types, ad apted from the four 



template types - E, Sbc, Scd, and Im - from IColeman et al.l (119801 ) . and the starburst 



templates SB2 and SB3 of iKinney et al.l (119961 ). were used. These templates are presented 
graphically in Figure [TJ These galaxy spectral types are derived from empirical data, and 
represent the range of non-AGN galaxy spectral energy distributions (SEDs) from redder 
to bluer colors; these starburst SEDs, for example, represent very blue galaxies not signifi- 
cantly reddened or obscured by dust. Interpolation was used between the six spectral types 
to produce a grid of 31 possible galaxy SED fits. 



Using the multiband photometry, the COSMOS photometric redshift code (IMobasher et al 



20071 ) was used to derive a photometric redshift z p for each galaxy, with 68% and 95% con- 
fidence intervals computed above and below that value. To determine the accuracy of the 
code, z p values were compared with spectroscopic redshifts (z s ) in 868 galaxies w ith z < 1.2 



and %ar < 24 with secure redshift measurements from the ^COSMOS survey (jLilly et al. 



20071 ). The galaxy spectral types (20% early-type, 63% spiral, 17% starburst) are evenly 
distributed with redshift in the spectroscopic sample, and about half of the sample (45%) is 
at z s < 0.5. 



The detailed statistics of the z p and z s comparisons are described in IMobasher et al. 



(120071 ). primaril y in terms of th e para meter Az = (z p — z s )/(l + z s ). As given in Table 4 
and Figure 5 of IMobasher et al.l (120071 ). the nominal dispersion between these photometric 
redshifts and spectroscopically determined redshifts was rms (Az) = 0.033 for non-AGN 
galaxies. With respect to the different spectral types, IMobasher et al. (l2007h showed that 
early-type, spiral, and starburst galaxies have rms (Az) = 0.034, 0.030, and 0.042 respec- 
tively. 

We chose to compute our luminosity functions using the l/-band data for this study. 
The saturation limit for bright objects in that image was V ~ 18.5, and the 3a faint- 
detection limit was V ~ 26.4; so we conservatively chose as o ur apparent magnitu de range 
19.0 < my < 25.0. Details of the photometry are described in ICapak et al.l (120071 ). A rest- 
frame absolute My magnitude, based on the derived z p and using the K-correction for the 
best-fit spectral type, was computed for each galaxy. Ground and space-based images of a 
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representative subsample of the faintest galaxies (My ~ —16 and fainter) in this sample are 
presented in Figure [2j 



3. Luminosity Functions 



As with all surveys which rely primarily on photometric rather than spectroscopic red- 
shifts, the application of the COSMOS galaxy sample to the derivation of galaxy LFs re- 
quires great care to account for both random and systematic errors in the redshift deter- 
minations. Just as methods of computing photometri c redshifts have evolved and improved 



(IKoolll986l : IConnollv et all 1199a iLiu fc Green Ill998l : Bemtez Il2000l : Mobasher et al.l 12004 



20071 ) , so too have techniques to quanti fy and compensate for the effects of relatively large 



20031 : iDahlen et al.ll2005l : IPerez-Gonzalez et all 120051 ). 



1996: 


Liu et al. 


1998: 


Chen et al. 



In this work, we adapt the method used in Liu et al. ( 19981), updating it with additional 



Chen et al. 



2003: Dahlen et al. 



components similar to those used in more recent studies (e.g 
20051 : iPerez- Gonzalez et al.l 120051 ) to reproduce the faint-e nd slope of the LFs of C OSMOS 
galaxies. Our strat egy is based on th e l/V max method (jSchmidt fc Green I Il986l ). which 
is well described by IChen et al.l (120031 ) as a maximum likelihood method for estimating a 
luminosity function without assuming any parametric form. We account for photometric 
redshift errors by weighting the galaxies as probability-smoothed luminosity distributions at 
the redshifts where they are measured. 



3.1. The modified 1/V max method 

Consider a galaxy with an apparent magnitude mj in a passband /, and redshift z ±<7. 
If a=0, then the absolute magnitude is 

M f = m f - 5 log(dr(z)) - 25.0 -k f (z) 

where di(z) is the luminosity distance in Megaparsecs, and k(z) is the K-correction at that 
redshift, in that passband, for the spectral energy distribution of the galaxy. The contribution 
of that galaxy to the luminosity distribution is then a delta function of amplitude unity at 
redshift z. 

In the case where a > 0, and the error distribution is Gaussian, the galaxy can be 
thought of as adding a series of fractional contributions to the luminosity distribution in the 
redshift space surrounding z. Such a fraction at, for example, redshift z + 5z and with a 
differential redshift width dz, would have an absolute magnitude 
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7 

and have an amplitude 



Ml = m/- 5 log(d L (z + Sz)) - 25.0 -k f (z + Sz) 



N z+Sz = P G (z + 8z, z, a)dz I A G (z + Sz, z, a) 

whe re Pr? and A G are the Gaussia n probability function and its integral, respectively (see, 
e.g. 



Bevington &: Robinson Ill992l ). 



For the photometric redshifts of the COSMOS survey, the redshif t error distribution 



is not Gaussian, but rather can be modeled as two half-Gaussians (ICapak et al.l 12007 



Mobasher et al.l 120071 ). where the 68% confidence interval on the lower and upper limits 
are ai and a u respectively. For a galaxy which a photometric redshift z p , the amplitude of 
the luminosity contribution function would be 

N 2+ fc = P G + Sz, z, ai)dz / A G (z + Sz, z, a t ) 



for z < z p , and 



for z > z p . 



N z+Sz = Pg(z + Sz, z, a u )dz / A G (z + Sz, z, a u 



This "fuzzing" of a galaxy's luminosity distribution in redshift space is straightforwardly 
achieved numerically, with a choice of dz « o to minimize random magnitude errors. For 
this COSMOS dataset, we used dzi = 0.02 ai and dz u = 0.02 a u . This divides each galaxy 
into a Gaussian-weighted luminosity distribution with 300 bins, from z—3ai to z+3a u . The 
entire distribution for each galaxy is normalized to unity. 

In the standard l/V max method, each galaxy contributes a weight to the luminosity 
function equal to the inverse of the accessible volume within which it can be observed. The 
accessible volume, referred to here as V max , is the total comoving volume within the redshift 
boundaries of the sample, where the given galaxy could be and fall within the selection 
criteria of the sample. In our case, the relevant criteria are the bright and faint apparent 
magnitude limits, and the effective solid angle of the COSMOS survey. 

In the case of a probability-weighted luminosity distribution for individual objects, it is 
straightforward to compute V max for each fractional galaxy; correspondingly, its contribution 
to the luminosity function is (l/V mtKC ) x N 2+< 5 2 . Assembling the luminosity function is then 
a matter of summing those contributions within absolute magnitude bins. 
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3.2. Redshift Limits And Sample Size 

Since the primary goal of this work was to examine the LF faint-end slope, the upper 
redshift boundary was determined mainly by our desire to sample with high completeness 
at least as faint as My ~ —16.5 in the entire redshift range. For a typical galaxy in the 
sample, depending on the galaxy's K-correction, this corresponds roughly to z < 0.4. But 
because each galaxy's luminosity is calculated as a probability-weighted distribution, there 
is a statistically significant contribution to the LF more than a full magnitude beyond the 
formal absolute magnitude limit. Thus, with the caveat that we are beyond that limit, 
we were also able to derive faint-end slopes of the LFs of galaxies in the redshift range 
0.4 < z < 0.5. 

Similarly, at the low redshift range, we are also able to measure fainter in absolute 
magnitude than the formally faintest detectable galaxy. However, to avoid large systematic 
magnitude errors and biases from structure in the local universe, we set a lower redshift 
bound of z > 0.02 for deriving the LFs. This meant that we have statistically meaningful 
luminosity contributions to the LFs down to an absolute magnitude of My ~ —12.8 for 
0.02 < z < 0.1. 

Within our apparent magnitude limits of 19.0 < V < 25.0, the COSMOS survey pho- 
tometric and z p catalog (ICapak et al.l 120071 ) contains 49,161 galaxies in the redshift range 



0.02 < z < 0.5. Below and above this redshift range, however, there are galaxies within the 
apparent magnitude limits whose probability-smoothed luminosity distributions contribute 
to the light within that range. Using 5 68 to denote the width of the 68% confidence interval 
for z p , we thus also include the contributions of all other galaxies whose luminosity distri- 
bution tails fall within 3 x 5 68 . (For example, a galaxy with z p = 0.51 and = 0.05 would 
contribute the portion of its luminosity distribution from 0.36 < z < 0.50, while a galaxy 
with Zp = 0.61 and 5qs = 0.05 would contribute from 0.46 < z < 0.50.) There are 31,659 
galaxies in the catalog that make such a partial contribution; even though many of those 
galaxies add only a tiny fraction of a galaxy into the z < 0.5 redshift range, we include them 
in our analysis for statistical completeness. Thus, a total of 80,820 galaxies are included in 
the sample used to derive the LFs. 

To check how the inclusion of these galaxies might affect the distribution of z v accuracy 
in the sample as a whole, we created a subset of the 80,820-galaxy sample where the galaxies 
have a <5 68 of at most 0.1. This subset of 41,237 galaxies (51% of the entire sample) effectively 
contains objects whose axe no more than 3 times the rms {Az) of the spectroscopically 
tested accuracy of the COSMOS z p code. We then computed the mean, median, and rms of 
the S 68 / (1 + z) values for the 41,237, 49,161, and 80,820 galaxy samples. We give the results 
in Table 1. As might be expected, the median, mean, and rms values all increase as the 
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constraints in redshift and are lifted and the sample sizes grow. The median Se&/(1 + z) 
value increases only modestly on an absolute numerical basis (0.042 to 0.043 to 0.052); the 
mean and rms increase more substantially. 



3.3. Simulations 

Computing LFs using "fuzzy" galaxies with photometric redshifts is clearly vulnerable 
to a set of systematic errors that would not be present for a galaxy sample with secure 
spectroscopic red shifts. In this work, we use the standard galaxy LF parameterization of 



Schechter I (119761 ). where M* is th e characteristic magnitude, <fi* the characteristic number 



density, and a the faint-end slope (ILin et al.lll996l ): 



0(M) = 0* x 0.41nl0exp (_i -°- 4 (m-a/,)) [ 10 -o.4(M-M,)(i+a)] dM 



(1) 



In any given redshift bin, objects near the peak of the LF — that is, near M* — will 
have part of their light distributed toward brighter magnitudes; and objects at the bright 
and faint end of the galaxy sample will have their light scattered still further. This would 
cause an overestimate of those parts of the LF which contribute the least light. Additionally, 
as the luminosity contributions of the galaxies are distributed across a large redshift range, 
light from galaxies inside a given redshift bin will sometimes be scattered out of that bin; 
and since lower- z redshift bins have smaller volumes than higher- z ones, there is the risk 
that more light would be added into lower-redshift bins than would be removed from them. 
This could bias the distribution of derived absolute magnitudes in those bins. 

We qua ntify and correct for these errors using Monte Carlo simulations in the manner 



described by iLiu et al.l (119981 ). For an arbitrary fixed M* and 0*, we created populations of 
galaxies that followed Schechter functions with faint-end slopes ranging from —0.2 < a < 
—2.2. For each a, we populated a simulated COSMOS survey volume with the corresponding 
galaxy population. We then "detected" these galaxies based on the magnitude and redshift 
limits of our survey; and for each galaxy detected, we randomly added an error to the redshift 
of that galaxy consistent with the measured dispersion of the COSMOS photometric redshift 
catalog, rms (Az) = 0.033. Finally, we assigned to that galaxy a 5 68 equal to the median 5 68 
of the 80,820-galaxy sample, i.e. 0.052 x (1 + z s ). The LF for this simulated galaxy sample 
was then computed using the modified l/V max method described above. 

We generated 150 such simulated LFs for each a tested. To illustrate the results of these 
simulations, we present three sets of them in Figure [31 The dashed lines show the input 
faint-end slopes (a = —0.7, —1.1, and —1.5) for the simulations. Of the 150 simulations, we 
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exclude the four outliers furthest above and four furthest below the input a; the remaining 
142 simulations (i.e. 95%) yield results that fall within the envelope bounded by the solid 
lines above and below each dashed line. The a represented by each of those solid and dashed 
lines is given in the figure. 

As Figure [3] shows, the "fuzzing" of the galaxies due to z p uncertainties does cause 
systematic biases of the calculated faint-end LF slope. At a ~ —11, the bias is negligible, 
and the fuzzing basically just increases the uncertainty of the measured slope. For steeper 
input a, however, the output a is quantitatively biased toward steeper values; and for 
shallower input a, the opposite is true. The bias increases as input values of a grow more 
extreme; an input of a = —1.9, for example, produces output a in a 95% envelope ranging 
from —2.3 < a < —1.8, whereas an input of a = —0.4 produces corresponding output of 
+0.1<a<-0.5. 

From all the simulation results, standard bootstrap methods were used to estimate the 
expected systematic offsets in a for an actually observed galaxy population. These offsets 
were then used to correct the computed luminosity functions, to recover the original faint- 
end slopes of the galaxy samples in each redshift bin. This correction affects the overall slope 
only, and is not intended to remove any inherent, non-power law structure in the observed 
LF. Also, this correction strategy was optimized to recover the faint-end slope, rather than 
the charact eristic magnitude M* or the number density normalization <ft* (as was done, for 



example, in lChen et al.ll2003l ). We thus do not attempt to measure those parameters in this 



work. Accurate measurement of those values for the COSMO S survey are presented , for a 



slightly different galaxy spectral type classification scheme, by IScarlata et al. I (120071 ) . 



3.4. Probability- Weighted Luminosity Functions 

The measured luminosity functions for the entire galaxy sample, corrected for this bias 
in a, are presented in Figure HI Only the faint ends of the LFs, operationally defined here 
as galaxies fainter than My ~ —18, are presented. LFs were calculated in five redshift 
bins: 0.02 < z < 0.1, 0.1 < z < 0.2, 0.2 < z < 0.3, 0.3 < z < 0.4, and 0.4 < z < 0.5. 
Additionally, we divided the galaxies into four subsamples according to galaxy spectral type: 
Type 1 (ellipticals), Type 2 (Sbc), Types 3 and 4 combined (Scd + Irr), and Types 5 and 6 
combined (low-extinction starbursts). For each subsample, LFs were also calculated in the 
same five redshift bins. The results are presented in Figure [5j 

There are three primary sources of errors in the LFs: (1) the systematic error in the LF 
slope and absolute magnitude determinations, described in the text above and characterized 
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using simulations; (2) a Poisson-like error which derives naturally from the modified 1/V max 
method, which is the reciprocal of the square root of the total number of fractional galaxies 
in each bin; and (3) a non-Gaussian error as a function of absolute magnitude, due to 
the asymmetric uncertainty in the photometric redshift determination of each galaxy. The 
second source of error, because of our large galaxy sample size, is much smaller than the third 
source of error. We computed the error from that third source with a standard bootstrap 
technique, using 150 random samplings (with duplication allowed) of the observed dataset 
to determine the 68% confidence intervals for $ in each magnitude bin. 

Each LF segment was fit to a power-law slope using weighted least squares. To avoid 
possible biasing of the faint-end slope by galaxies near M*, which as we mentioned above 
is not well determined by our technique, we only use data fainter than My = —18.3 in our 
fits. In one case, we make a more restrictive magnitude cutoff: the lowest redshift bin in the 
starburst spectral type. There, the measured LF drops discontinuously at My = —17. We 
suspect this has occurred because M* may be quite faint for this galaxy type and redshift 
bin, thus distorting our measurement of the LF brightward of that point. So to ensure 
that the a fit to that LF is not correspondingly distorted, we use only data fainter than 
My = —16.8 for that particular fit. In both Figures @] and El the measured LF data points 
are plotted as symbols, and the best-fit a values are plotted as dotted lines. The faint-end 
slope fits are summarized in Table 2. 



4. Discussion 



Because of the combined large area and depth of the COSMOS survey, the luminosity 
functions presented in this work provide a glimpse of a across a substantial range of redshift 
in a single, consistent dataset. With this view, our results show that for all spectral types 
combined, a = 1.24 ± 0.07 for the local (0.02 < z < 0.1) universe. As the redshift increases, 
a flattens out somewhat, and is —1.12 ± 0.10 in our highest redshift bin (0.4 < z < 0.5). 

Our local LF is consistent with results from the two largest local galaxy surveys to date, 
which have comparable (~ 10 5 galaxies) sample sizes to our study here. The 2 dFGRS survey 
found that, for the bj band galaxy luminosity function, a = —1.21 ± 0.03 (jNorberg et al. 



2002h and -1.18±0.02 for the redshift range 0.02 < z < 0.25 JCroton et al.lboosh . The SPSS 
has a = — 1.05 ± 0.01 for a slig htly redshifted r-band galaxy LF for g alaxies brighter tha n 
M r ^ — 17 (IBlanton et al.l 120031 ). and a ~ —1.3 at fainter magnitudes (IBlanton et al.l 120051 ). 



Blanton et al.l (120051 ) has further shown that, with the appropriate conversion of the 2dFGRS 
bj data, it and the SDSS g-band LFs have consistent low-l uminosity slopes. At higher 



redshift, our results are consistent with IScarlata et al. I (120071 ). who independently derived 
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oi = —1.26 ± 0.15 in the range 0.2 < z < 0.4 for a portion of the COSMOS survey area. Our 
results are also consistent with results from oth er surveys given in the literature - for example, 
with the V-band LF derived from the VVDS Jllbert et al.ll2005h . where a = 1.21 ± 0.04 in 
the range 0.2 < z < 0.4. 

The formal errors in our a measurements are higher than those of most of these other 
studies; this is probably mainly because of the systematic slope uncertainties that we attempt 
to account for with our simulations. Our overall agreement, however, appears to confirm 
that we have properly accounted for the errors that result from representing galaxies as 
probability-smoothed luminosity distributions. 



4.1. LFs as a Function of Galaxy Spectral Type 



The luminosity functions in the four galaxy spectral type bins we used - SED templates 
of early-type, Sbc, Scd+Irr, and low-extinction starbursts - follow the well-known pattern of 
steeper a for bluer galaxies. In our low-redshift bin, a increases from —1.10 ± 0.08 in early- 
types to —1.88 ± 0.18 in starbursts; the trend continues with increasing redshift, showing a 
similar steepening of a, from —0.52 ± 0.20 to —1.27 ± 0.15. 

As with the full galaxy sample, these type-specific results are consistent with the find- 
ings of previous work in the literature. Of course, due to differing galaxy selection criteria 
and redshift binning, exact comparisons are not always possible. Generally speaking, how- 
ever, for local galaxies, our red/early-type and intermediate spiral galaxy LFs are consistent 
with SDSS and 2dFGRS results; and previous samples of z < 0.1 galaxy populations have 
also shown very steep a for the bluest a nd most irregular galaxies: a = —1.87 for the 
CfA Redshift Survey (IMarzke et al.l 119941 ) , —1.84 for the Las C ampanas Redshift Survey 



(IBromlev et al.lll998T). -1.81 for th e SSRS2 (IMarzke et a l. 1998 ), and -1.9 from the Deep 



Multicolor Survey ( Liu et al. 1998h and SPSS (jNakamura et al 



comparisons with COM BO-17 (IWolf et al. 



2003) At h igher redshifts, 
20031. VVDS JZucca et aHl2006h . and with the 



COSMOS survey itself (jScarlata et al. 1 120071 ) show broad consistency across the various 
galaxy type and redshift intervals. 



4.2. a vs. z: Evolution Or Selection? 

In the context of the broad consistency of our results with those in the literature, perhaps 
the most striking result in this work is the clear trend, with every galaxy spectral type, of 
a flattening of the faint-end slope with increasing redshift. From our lowest redshift bin to 
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our highest - i.e., from z ~ to z ~ 0.5 - Aa = 0.58,0.24,0.35, and 0.61 respectively for 
early-types, Sbc, Scd+Irr, and low-extinction starbursts. 

On the surface, this trend may not appear to be consistent with previous work. Much 
of the work to derive the evolution of galaxy LF parameters M* and $* as a functio n 



of redshift (e.g., iLin et al.lll999l : IWoff et al.ll2003l : Baldrv et al.ll2005l : IWillmer et all 120061 ) 



in fact, depends upon the assumption that a does not evolve with redshift, or at most 
weakly evolves to z < 1. We have focused here on the measurement of a rather than those 
other parameters; that, and the fact that all of our measurements have come from a single 
dataset, plus the fact that each individual determination of a is consistent with previous 
work, supports the likelihood that this observed flattening trend is real. 

The question is, do these changing slopes represent true evolution in a, or do they 
reflect our ability to detect different galaxy populations as a function of redshift? The l atter 
possibility can be dis cussed in the context of, among others, Ide Lapparent et al.l (120031 ) and 
Blanton et al.l (120051 ). who suggest that the faint end of the field galaxy LF is comprised of 
a composite population of dwarf and non-dwarf galaxies, each with its own functional form. 
This would mean that a s i ngle p ower law is not quite sufficient to describe the LF faint end 
accurately. iBlanton et al.l (120051 ) further suggests, through detailed examination of the faint 
galaxy population in the SDSS, that a large fraction of these dwarf galaxies may have very 
low surface brightness, and are thus not included in most faint end LF measurements. 

Our dependence on photometric redshifts places an important caveat on the interpre- 
tation of our data: by using "fuzzy" galaxies, any second-order deviations from a power law 
at the faint end of the LF may well have been smoothed out, and are thus not recoverable 
from our LF measurements. So if a deviation from a single faint-end power law does occur 
at very low luminosities, we cannot address that issue with this work. 

Due to the substantial depth of the COSMOS survey imaging, it is likely that we have 
successfully measured a larger fraction of low surface brightness dwarf galaxies than wider- 
area, shallower surveys such as SDSS or 2dFGRS. The steepness of our low-redshift LFs 
may reflect this. Even the COSMOS survey depth, however, does not allow us to measure 
a fainter than My ~ —17 at z — 0.5, so we cannot say if the flattening of a in our higher 
redshift bins is due to the non-detection of these dwarfs. There may be some circumstantial 
evidence, however, to support that picture. For example, iDahlen et al.l (120051 ) measured for 
the GOODS survey a = —1.37 for the rest-frame B-band from 0.1 < z < 0.5, which is 
somewhat steeper than most LF measurements in this range. However, the GOODS survey 
is very deep; so in this broad redshift bin, a large detected fraction of faint, low surface 
brightness dwarfs near z > 0.1 could have driven a to a steeper value for the full range. 
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As the spectroscopic portion of the COSMOS survey (ILilly et al.l 120071 ) is completed, 
we will be able to address this question in more detail, as we deconvolve the faint end of the 
galaxy LF as a multivariate function of color, morphology, luminosity, and redshift. 



5. Conclusions 



Using the COSMOS multiband photometry and photometric redshift catalog, we have 
constructed faint-end rest-frame V-band luminosity functions for the galaxy population at 
0.02 < z < 0.5 in the COSMOS survey volume. Since we are using photometric redshifts, 
we have computed these LFs by treating galaxies as weighted probability-smoothed lumi- 
nosity distributions, and using a modified l/V max method. A total of 49,161 galaxies have 
photometric redshifts that fall in this redshift range; within and outside this range, a total 
of 80,820 galaxies contribute to the derived LFs. Extensive Monte Carlo simulations were 
used to characterize and account for the systematic and random errors of this technique. 

For all galaxy spectral types, the LF slope ranges from —1.24 to —1.12 from the lowest 
redshift bin to the highest. In the lowest redshift bin (0.02 < z < 0.1), where the magnitude 
limit is My < —13, the slope ranges from a = — 1.10 ± 0.08 for galaxies with early-type 
spectral energy distributions (SEDs), to a = —1.88 ± 0.18 for galaxies with low-extinction 
starburst SEDs. In each galaxy SED category (early- type, Sbc, Scd+Irr, and starburst), the 
faint-end slopes grow shallower with increasing redshift; in the highest redshift bin (0.4 < 
z < 0.5), a = —0.52 ± 0.20 and —1.27 ± 0.15 for early-types and starbursts respectively. 

All of our derived type-specific LFs, across our redshift ranges, are broadly consistent 
with the findings of previous authors. Our results thus show a flattening trend for a with 
increasing redshift for each spectral type. It is unclear, however, if this is evidence of a 
evolution in the galaxy LF, or of preferential selection of dwarf galaxies in the local universe. 
The steepness of a at lower redshift could be qualitatively explained, for example, by large 
numbers of faint dwarfs, perhaps of low surface brightness, which are not detected at higher 
redshifts. We will address this question in a future paper, as the full set of COSMOS data - 
in particular, spectroscopic redshifts - are obtained. 
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Table 1. Galaxy samples and z p confidence levels 



# Galaxies 


Redshift constraint 


(568 constraint 10 


median 568 /(I + z) 


mean <5e8 /(I + 2 ) 


rms 568/(1 + z) 


41237 




<5 6 8 < 0.1 


0.042 


0.041 


0.053 


49161 


0.02 < 2 < 0.50 




0.043 


0.050 


0.073 


80280 






0.052 


0.081 


0.239 



a Galaxy spectral types as defined by Mobasher et al. 
b <568 = the width of the 68% confidence interval for z p . 
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Table 2. Luminosity Function Slope Fits 



Galaxy Spectral Type a Redshift Range Slope a 



T1-T6 (All) 


0.02 < z < 0.1 


-1.24 ±0.07 




0.1 < z < 0.2 


-1.18 ±0.07 




0.2 < z < 0.3 


-1.09 ±0.08 




0.3 < z < 0.4 


-1.12 ±0.08 




0.4 < z < 0.5 


-1.12 ±0.10 


Tl (Ell) 


0.02 < z < 0.1 


-1.10 ±0.08 




0.1 < z < 0.2 


-0.81 ±0.09 




0.2 < z < 0.3 


-0.60 ±0.12 




0.3 < z < 0.4 


-0.53 ±0.16 




0.4 < z < 0.5 


-0.52 ±0.20 


T2 (Sbc) 


0.02 < z < 0.1 


-1.16 ±0.07 




0.1 < z < 0.2 


-1.13 ±0.07 




0.2 < z < 0.3 


-1.01 ±0.10 




0.3 < z < 0.4 


-0.75 ±0.15 




0.4 < z < 0.5 


-0.92 ±0.16 


T3+T4 (Scd+Irr) 


0.02 < z < 0.1 


-1.46 ±0.07 




0.1 < z < 0.2 


-1.37 ±0.08 




0.2 < z < 0.3 


-1.28 ±0.09 




0.3 < z < 0.4 


-1.19 ±0.10 




0.4 < z < 0.5 


-1.11 ±0.15 


T5+T6 (Starbursts) 


0.02 < z < 0.1 


-1.88 ±0.18 




0.1 < z < 0.2 


-1.65 ±0.14 




0.2 < z < 0.3 


-1.53 ±0.10 




0.3 < z < 0.4 


-1.35 ±0.11 




0.4 < z < 0.5 


-1.27 ±0.15 



a Galaxy spectral types as defined bv lMobasher et al.l ll2007i ) 
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Fig. 1. — Spectral energy distributions from iMobasher et al 
MOS galaxy spectral types and photometric redshifts. 



(120071 ) used to compute COS- 



-21 - 



22M4HI iik= 2J.II r= ilJU M- IUC7 



t 





« 



iu-iJ.17 *=(ij47 M=-LLiO m= U.7U *= d.lH M= -UJM 




# 1 \ I 

^>:oj m-^j.r,i i-JLiu M— u.'w 



22HSflI4 m=24.7J i=4Ui M=-]}.HJi 



■ 3iV 




icuim iu=ujii »=imk? l^-.isjr? 

i • / 

* J 

2:51343 m= 24.411* i=O.HJ \t= JJ.Jfe 




179334)13 *4S? M= 14.17 




if * 

2792! >S m= 244)4 j= U.t 1 11= , I5ii7 



mwsmtmmmwmmmn^mwn mm. J* w < 



KL237H HI-I.V7J j=.IJ.t? M=-];.7J 



■ 



lriSfriv7> ni=llj;.l nll.lJ* M=-NJ4 




I DO.) 1 40 ni-i?ji|t i-L'Jii.\l- H.JM 




l>IW*7II dkMJI IMS 11= .IftJlt 




IWSH3H m=>2.55 1=8 4N Vl= .Ifi.U IV29972 ni=23»H I=\t2i 3]=.1<1JJ 



Fig. 2.— Hubble ACS J 8 i 4 and Subaru B images of a representative selection of very 
low-luminosity galaxies in the COSMOS field. Each image is 15" across, and each image 
pair is labeled by the COSMOS catalog ID number, photometric redshift, and apparent and 
absolute V magnitudes of the object at the image center. 
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Fig. 3. — Examples of the results of simulations testing the biasing effect of probability- 
weighted luminosity functions on the measurement of a. The dotted lines represent a of 
the input LFs, and the solid lines show the boundaries wherein 95% of the output LFs are 
contained. The three numbers above each set of lines are the a of the upper bound, input, 
and lower bound LFs respectively. 
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Fig. 4. — Faint-end portions of the V-band galaxy luminosity functions, $(My), in the 
redshift range 0.02 < z < 0.5. Each LF is offset by a constant for clarity. The best weighted 
least-squares-fitting faint-end power law slope for each LF has been overplotted (dotted lines) 
and labeled with its redshift bin and slope. 
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Fig. 5. — Faint-end portions of V-band luminosity functions, $(My), for galaxies divided 
by spectral type. Symbols are the same as in Figure HI 



